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ABSTRACT 


‘An acoustic ray tracing algorithm is developed and coupled with a 
thermodynamic upper ocean mixed laver model. For a test case, the coupled mixed 
laver-acoustic model is applied to a specific area in the western Mediterranean Seca. 
Climatological atmospheric forcing is used to provide boundary conditions for the 
nuxed laver for short periods of time (from few hours to three days) during different 
seasons. The response of the acoustic model to the predicted changes in the sound- 
speed profile is analyzed to show dependence of acoustic propagation upon the surface 
atmospheric forcing and the season. The atmospheric factors such as wind, rain, and 
solar irradiation have almost no effect on the propagation of rays emanating from a 
deep transmutter. In the case of a shallow source, the wind is the most dominating 
factor wnich influences the acoustic propagation. The effect of heavy rain with light 


Mind 1s also examined. 
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lL. iNTRODUCTION 


An acoustic rav tracing program 1s coupled with the Oceanic Boundary Laver 
Model developed by Garwood [Ref. 1]. The OBL model is a one-dimensional. second- 
order turbulence closure, vertically integrated model of the ocean surface turbulent 
boundary layer, using a two-component turbulent kinetic energy budget with a mean 
turbulent field closure. 

Fisher [Ref. 2] investigated the variability and sensitivity of a coupled model 
svstem. He found that the OBL model, when integrated in time at a single point 
(Ocean Station Papa 50°N 145°W), predicted mixed-layer structure better than did the 
Expanded Ocean Thermal Structure (EOTS) svstem which was currently in Usem@esrme 
Fleet Numerical Oceanography Center (FNOC). McManus [{Ref. 3] evaluated the 
acoustic performance of a coupled model system at a line of stations in the northeast 
Pacific Ocean. In both cases, the thermodvnamic model was initialized with observed 
temperature profiles, and the surface boundary conditions were given by the currently 
available metcorolegical informations. Then, the thermodynamic forecasts were input 
Into acoustic models, such as RAYMODE or FACT, and the acoustic performance 
Was analyzed using the median detection range (MDR) and the convergence Zone 
range { G7 

This research 1s the first attempt to link in a single program the OBL model with 
an acoustic model. Simplicity (compared to the operational models available in the 
U.S. Navy) and classification restrictions lead us to develop in Chapter II an algorithm 
for acoustic ray tracing. As no such routine was available at NPS, a copv is attached 
In the appendix for further use by students of the Air-Ocean Sciences Department. This 
simple subroutine allows the mfluence of the atmospheric forcing on the underwater 
sound propagation to be qualitatively analyzed. A summary of the leading principles 
and equations of the OBL model ts given in Chapter III. Chapter IV develops in detail 
the actual coupling of the two models into a single computer code. Chapter V gives an 
exampie of the use of this coupled model applied to a specific area in the)“ meme 
Mediterranean Sea, for different periods of the year having particular acoustic 
propertics. All the stmulations were integrated in time using climatological data over 


short time periods varving from a few hours to three days. 
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If. ACOUSTIC MODEL 


A. INTRODUCTION 

Since we assumed in this research that the ocean 1s horizontally stratified. the 
temperature I and the salinity S are only functions of the depth y and cannot vary 
with the range z. Hence, the speed of sound c 1s onlv a function of depth y. In that 
Siocwaccorcing to Ziomek {Ret. 4: p.236], the general form of the equation for the 


horizontal range travelled bv an acoustic ray is the following : 


v cc 
= —_ cr oS en d< us 
i \, (eri ey 


where b is the ray parameter and is given by 
b = sinB(yq) / c(¥9), aD 


(¥9.Z)) are the coordinates of the source, and B(¥y) is the initial angle of propagation. 
Thus, theoretically, knowing the sound-speed profile versus the depth, we can 
plot the curve giving ithe path of an acoustic ray. 
We chose not to use equation (2.1) for the following reasons: 


According to Snell’s law, 


sinB(y) _ sinB(¥y) _ b (2.3) 
c(¥) e(¥p) - 


and, at a turning point. BOvtp) = 7/2 and CY tp) = 1/b such that the denominator of 
the integrand of (2.1) goes to zero and the integration cannot be carried out. Thus 
(2.1) 1s only valid between two turning points. 

Numerically. the integration of (2.1) would be carried out between all the turning 
points of the ray by discretizing the sound-speed profile c(v) with as small a depth 
increment as required to get an acceptable result. Furthermore, during the integration, 


we would have to test for the occurence of a turning point. 


te 


However, a sound-speed profile can be approximated by straight line segments 
matching the profile as best as possible. The more segments chosen, the better the 
computation. Roughly, we often choose 10 to 15 segments depending on the shape of 
the profile. Then, as explained in the next section, the integration of (2.1) can be 


solved analytically for each segment of straight line and leads to very simple equations. 


Ze 
' 
! 
I 
; 
' 
1 
i 
| 
I 





Figure 2.1 Ray path confined to the YZ plane. 


We shall keep in mind that this acoustic model has to be coupled with an 
Oceanic Planetary Boundary Layer (OPBL) model [Ref. 1] which resolves the 
temperature, salinity, and depth of the upper mixed layer of the ocean using a one- 
meter-step increment for the depth. 

These are the reasons why we chose to discretize the sound-speed profile with a 
one-meter-depth increment, assume the profile to be a straight line segment within each 
depth increment, and we use the equations derived in the next section instead of 
numericallv evaluating the integral of (2.1). We found this method easier for handling 


the turning point problem. 


B. SPEED OF SOUND AS A FUNCTION OF DEPTH WITH CONSTANT 
GRADIENT 


The sound-speed profile 1s given by 


CO) a= Co + BX (Y-¥p) (2.4) 
12 


I 
where g 1s a constant (with units of s” ) referred to as the gradient since 





dc(v) 
es aS 
dy (2.5) 
Starting with Snell’s law (2.3), differentiation leads to: 
cosB(y) dB = bdc(y) = be dy. (2.6) 
Referring to Figure 2.2, it can be easily seen that: 
= cosf(y) (267) 
re y). . 
P da a 
dy 
ye 
Figure 2.2 An infinitesimal element of arc length ds 
at an arbitrary point P along a ray path in YZ plane. 
Substituting (2.7) in (2.6) vields 
AB(y) 
~= bg. 8 
ip bg (2.8) 


IS 


Equation (2.8) indicates that the curvature along a ray path is constant. Thus. the rav 


path is an arc of a circle. From Figure 2.2, we have 


dy 


dz = tanB(v) dv = sinp(y) Sule 


Using (2.6), we get 


ESD 


o 


dz = sinB(v) 


Integration leads to 





(2.9) 


(2.10) 


(2.11) 


where B, = B(y,) is the initial angle of propagation at the source, and § = fi(y) is the 


angle of propagation at position P. This leads to 





l 
= Ors (c(v) - Co), 


and, using Snell's law (2.3), we get 





: 
¥ = Yo + > (sinB - sinBy). 


(2.12) 


(2nies 


(2.14) 


We can verify that (2.12) and (2.14) are the parametric equations (parameter f) of the 


Cliclercivenvine: 
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centered at { Vp - sinBy. Dag 4s cosBy bg) ) and having a radius of 1 bg. 


C. ALGORITHM 

For this research, we used a modified version of the basic FORTRAN subroutine 
RAY given in appendix A. In this section we are going to analyze in detail the 
different parts of the algorithm used in the subroutine RAY. 

The goal of this subroutine is to plot the acoustic ray path, that is, the curve 
giving the range z in kilometers versus the depth y tn meters of an acoustic ray 
emanating from a source at a depth yp at a given initial angle of propagation. The ray 
path plot is presented beside the graph of the sound-speed profile c(v) in misec versus 
depth y in meters. 

As we mentioned previously, we used a one-meter-step depth increment and kept 
track of the depth all along the ray bv using an increasing or decreasing index k 
depending on the rav going upward or downward. The index k is, in fact. the integer 
value of the depth. 

The angles of propagation B and Bo are releremecd from te Y axis, For plotting 
purposes, we defined the depth vector { y,, 
me Ny}. At the th point of the rav path. the depth index k ts 


1 = O}XN >} and the range vector { ae 1 = 


the “initial” angle of propagation is By. the sound-speed is c, and, its gradient. gy), is 


given by 


> 


2k ~ “k+1~ °k: Gale 


Avoine (1+ jth point on the ray path. the depth index is k£I1 the “final” angle of 
propagation is J, the sound-speed is c, +; and the gradient is g, 4.) . An illustration 


of these different parameters is given in Figure 2.3. 
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For that case, according to (2.3) and (2.12), we would have 


eh nem! (2.18) 

and 
= >OS 2.19 
Zi+] = Zi as om (cosB, = cos) ( ) ) 
=K-] 

with 

sinh  sinf 

Chall “source 


O (Grae Ck Cri 





Figure 2.3 Parameters describing a one-meter-depth increment 
positive gradient sound-speed profile 
for an upward travelling acoustic ray path. 
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Now we are going to examine in detail the different cases occuring in the 
Subroutine neading to some IF... THEN...ELSE... statements, depending on the sign of 
the gradient (or the curvature) and the upward or downward direction of the ray at the 
jth point. 

l. Case 1: g, ,>0, ¢, > 0. 
At the ith point of the ray path P,, corresponding to the depth index k, the 


gradient above is (see Figure 2.4) 
Bees Ck Sk] (2.21) 
and the gradient below 1s 


mec ok Oe 


Cw, on Cee C 





Figure 2.4 Case |: 2, _,>0, 2, > 0. 


We have three different cases to consider depending on the value of B, with 
regard to m/2 and B.. The critical angle B. is the value of the initial angle of 
propagation at P, leading to a final angle S = m/2 one meter deeper, that is, leacing to 


a turning point one meter deeper. 
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According to Snell's law, 


| Ch : 
Se)b es — (2255 
Sheed 
Or 
Se Re | 
B. = arcsin ——. (2.24) 
k+1 


Let us examine the three cases B) > 2/2, B) <B, and B. <P) < 7/2 since a strict 
equality does not apply for real numbers in FORTRAN language. 
Cee 2 
According to (2.3) and since the ARCSIN function gives a result between 0 


and 7/2, the angle of propagation at P;, , is: 
J = m - arcsin(bc,_7) - (22255 


Substituing for B in (2.12), we get the coordinates of the next point Bea 


Vig y = kel (2.26) 
and 
| 
eR ao be, oho + coslarcsin(be,_1)]} (222) 
with 
b = sinB source ! Ssource - (2.28) 


5. By < Bo: 


This time, using the same equations (2.3) and (2.12) at the point P; 4). we 





have : 
SB = arcsin(be, 4 4), (2529) 
and 
eine in(b one 
Ai ee oe be, (cos, - cos[arcsin(be, 4 1)}} . (37) 
z oir 
oc Le < By = Toes 
in this case, we have to deal with a turning point within the segment 
ED. 
oe 


Let us derive the general formula for computing z: 4, When a turning point 
ieeencountered. [his formula will be used when we study the cases generated by the 
other sign possibilities of the gradient g. 


Applying (2.12) first from P, to P’. then from P’ to P we have (see Figure 





Eo) : 
, i | : 
We Ly sas (Cosi, cosn; 2) (232) 
and 
i, i be [cosm'2 - cos(m - By], 238) 
or 
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I 
zZ= 2) + ae (2.34) 
and 
| 
Z = 2 eo 535° 


bg 


Substituing (2.34) in (2.35) yields 


Z = Zo + Seok ae (2.36) 






\ meter 


Figure 2.5 Turning point treatment. 


Going back to our case of interest and using (2.36), we get the coordinates 


Vie =y=k (2.37) 


and 





2 ; 
7 el AR — cosy . (33) 


faCase 2 : Eke] < (), g). <0 
This case is very similar to case 1, but now the critical angle f. is given by : 


pe arcsin(c,/Cy_1) . (2739) 


Following the notation of Figure 2.6 and applying the same basic equations (2.3), 
t2)!2), and (2.36), we have three new cases to exanune. 





Figure 2.6 Case 2: 9, 1 <0, 2, <9. 


G@ Py < mi2: 


b= aresin(be, + 1), (2.40) 


al 


= hoe (2.41) 





and 
I 
Za quae bg; (cosP, - cos{arcsin(bc, 4 1)}} . (2.42) 
ee ee 
B= 1 - arcsin( be, 1) , (2.43) 
Vee ole (2.44) 
and 
: = 7. + i f sB “s in(b BAS 
Zea Z; ba, je g + coslarcsin(be,_,)j} . (2.45) 
CU eae 
Se (2.46) 
and 
: 2 ; 
Zi+] — Z + be, oro : (2.47) 


Coa 


9 > oO, < 
Case 3 oe 
According to Figure 2.7, we now have only two separate cases to consider, 
depending on the value of B, with regard to m,2. Returning to the two previous cases, 


we will find some similarities. 


tl 
ty 


ds Poe 7/2: 
This case is similar to case 1.a, and equations (2.25), (2.26), and (2.27) 
apply. 
b. By < n/2: 
This case 1s similar to case 2.a, and equations (2.40), (2.41), and (2.42) 


apply. 


' 
! 
| 
| 
| 
! 
1 
' 





PigUiee) ee ASeao ey Oey — 0, 


4. Case 4: g,_) <0, 2, >90 
According to Figure 2.8, we now have four different cases to examine with 


two possibilities for Pi... 


aresin —E (2.48) 
k-] 
or 
Ve arcsin —K (2.49) 
Keto 


2S 


The four different cases shown above reduce to cases previously analyzed in 


sections for case | and case 2. 


Yee 


T >B,> By 


Sear? het 





Figure 2.8 Case 4 open tney 20. 


we tba 
This case 1s similar to case 2.b, and equations (2.43), (2.44), and (2.45) 
apply. 
be eee tc 
This case 1s similar to case 2.c, and equations (2.46) and (2.47) apply. 
c. By < Be: 
This case is similar to case 1.b, and equations (2.29), (2.30), andaizaees 
apply. 


d. 1/2 > By > oe 
This case corresponds to case lI.c, and equations (2.37) and (2.38) apply. 
This last case ends the cascade of IF...THEN...ELSE... statements that we used to 


Write this subroutine. 


24 


5. Surface and bottom reflections 
The last point we need to discuss 1s how to handle surface and bottom 
reflections, that is, when the depth index k reaches the values of 0 or N. The value N 
corresponds to the maximum depth given in the profile and 1s assumed to be the depth 
of the sea floor. In the case where the sound-speed profile does not extend to the 
bottom, it is always possible to extrapolate the sound-speed profile to the bottom using 
the gradient of the profile corresponding to the last straight line segment or, using an 


average value given by the climatology. 


See cases |.2and 1.3 


See case 24 





Figure 2.9 Perfect surface reflection possibilities. 


a. Surface reflection 
In this simple acoustic model, we assume a perfect surface reflection. Thus, 
when the depth index k reaches the value 0, We just have to maintain symmetry (1.e., 
the angle of reflection is equal to the angle of incidence) with regard to the horizontal 
to get the correct angle By before applying the same equations as in cases 1.b, l.c and 
eeevol Figure 2.9 . 


a> 


b. Bottom reflection 
Following the same idea, we assume a perfect bottom reflection. Thus 
when the depth index k reaches its maximum value N, we maintain symmetry (1.e., the 
angle of reflection is equal to the angle of incidence) with regard to the horizontal (see 
lietire Zen): 
Notice that the fictitious gradients g_, and gx have been created in the 
program to handle the cascade of IF... THEN...ELSE... statements. 
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Figure 2.10 Perfect bottom reflection possibilities. 


D. SUMMARY AND DIRECTIVES TO USE THE SUBROUTINE RAY 
1. Main characteristics 
The computation in the subroutine RAY is based on representing an arbitrary 
sound-speed profile by straight line segments in one-meter depth increments and thus 
allows the use of any complicated sound-speed profile to obtain a precise and smooth 


acoustic ray trace. 
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A plot of the sound-speed profile is provided at the same time as the plot of 
ine Oundle of rays, 
The main assumptions underlying this program are the following : 
® perfect surface reflection, 
© perfect bottom reflection, 
¢ flat bottom. 
2. Features of this subroutine 
This subroutine uses many DISSPLA graphics statements. Hence, the main 
program calling subroutine RAY has to be executed using the command DISSPLA. A 
choice of multiple display devices is provided to the user through the use of comment 
fords im front of the statements CALL COMPRS, CALL TEK6I8, or CALL 
CX31(4107). 
The variables required when using RAY are the following : 
*) 40: depth of the source in mieters. 
pe) aumber of ray path desired, 


e BET: array (size M) of the initial angles of propagation B in degrces, 


source 
e BO: imtial angle of the upper ray of the bundle, 
See OB increment of muitial angles in BET, 
® MM-+1: number of points in the provided sound-speed profile. 
* CC: provided values of sound speed in m‘sec, 
Sey < Corresponding depth in meters, 
* RANGE: maximum range desired in kilometers, 
oe SN index of range, 
Sn otra. (Size NaN) Of Cepth im meters: 
® 7Z: array (size NN) of horizontal distances in kilometers, 
e WN: integer value of the depth of the sea floor in meters, 
Sm C : array (size N) of sound spced (m’sec), 
Seer oiray (Size iN} ofits cradient (sec”!), 
ome © = array (size \) of deptns (meters). 

The bundle of rays to be traced is defined by the number of rays M, the initial 
angle BO of the upper ray of the bundle , and the angle increment DB between two 
rays. inside the subroutine, we can aiso set different equations to define a bundle of 
rays. And finally, bv deleting these few lines involved with these computations, we can 


provide our own array BET of initial angles Boyce: 


No relation has been derived between NN and the maximum range RANGE 


since this relation depends on many different parameters such as the B ‘S.liee 


source 
sound-speed profile, and the maximum range. If during a plot a ray ends before 
reaching the right side of the graph (especially the steepest rays), a higher number NN 
has to oe set. For instance, a value NN =5000 seems to nandle most of the reasonable 
steep B.ource given for a4 maxinium range of 30 kim. 

The arrays Y and Z have been created for computing and plotting the ray 
paths and the arrays C, G and YC are for plotting the sound-speed profile. 

3. Example 
A result of the RAY subroutine is displayed in Figure 2.11. The provided 


parameters and sound-speed profile are listed in Table 1 . 














TABLE 1 


PARAMETERS AND DATASCSED (ie be ee 
OF SUBROUTINE RAY SHOWA rite UR Gel 


x0 ou: oe 0. ce 1482. 

M 10 30, 1352.9 
NN 5000 Be): 1484.1 
RANGE 30. 90. 1480.5 
N 300 hs) 1482.0 
MM 6 oS. 1484.1 


1486.2 


(WH) JINYY 





GO Wo Zari 0° O8yl 


(S/W) dSS 


(W) HId3d 


Figure 2.11 Ray tracing provided by subroutine RAY. 


I. OCEAN MINED LAYER MODEL 


A. INTRODUCTION 

The study of the oceanic turbulent boundary layer is a relatively recent field in 
Physical Oceanography. The model used for this research has been developed by 
Garwood [Refs. 1,5]. These papers give the formulation of a unified mathematical 
model of the one-dimensional, non-stationary, oceanic turbulent boundary layer. 

The study of these top few tens of meters of the ocean is Of C@msidenige 
scientific interest. It influences and can be related to the general circulation of the 
ocean [Ref. 6]. The thermal structure associated with this boundary layer must be 
considered when making medium and long-range weather forecast, since a large part of 
the atmospheric energy supply comes from the heat exchanged with the ocean. This 
laver is also a region of primary biological productivitv, which is of significant 
ecological and economic importance. Finally, as an important military application, this 
study can be used tn the modeling of acoustic propagation in the ocean, which is the 
goal of this research. 

1. Characteristics of the oceanic mixed layer 

The oceanic mixed laver is that fully turbulent region of the upper ocean that 
is bounded above bv the air-sea interface and below by a dynamicalv stable water 
mass. The wind and intermittent upward surface buoyancy flux through the surface 
(surface cooling at night or in Winter time for example) are the sources of mechanical 
cnergy for the generation of these turbulences. Figure 3.1 gives a general picture 
idealizing density and mean velocity profiles of the ocean mixed layer.! 

-. Generalities on the dynamics of the mixing 

The depth of the ocean wind-mixed surface layer is typically on the order of 
ten to one hundred meters. The horizontal scale size is that of the internal Rossby 
radius, typically 20 to 50 km. These two dominant scale sizes are much smaller tnan 
the horizontal scale size of the driving meteorological disturbances, water mass 


features, and distances to lateral boundaries. The approximation of local horizontal 


'In this chapter, in order to be consistent with the notation used in [Refs. 1.5] as 
in most of the geophysical sciences publications, we chose z to be the upward vertical 
axis and (x,y) to be the horizontai coordinates as depicted in Figure 3.1. The vertical 
component of the fluid velocity will be w and the horizontal components will be u and 
ne 


homogeneity for all mean variables is usually accurate and is a basic assumption in this 


model. 





Figure 3.1 Idealized model for ocean mixed laver.. 


A sharp density discontinuity of thickness 6 (see Figure 3.1) separates the 
layer from a stable non turbulent thermocline. Minimal stress at the bottom together 
with high turbulence intensity leads also to an approximate Vertical uniformity in mean 
velocity and density. We shall note that only small gradients in these mean variables 
give rise to large turbulent fluxes. 

The mechanical energy budget for the ocean mixed laver is depicted in Figure 
3.2. Deepening of the mixed laver is accomplished by entrainment of the more dense 
underlying water into the turbulent region above. This process leads to a potential 
energy increase and cannot take place without an energy source: the turbulent kinetic 
energy of the mixed layer above. 

Retreat occurs when the vertical component of the turbulence 1s insufficient to 
transport heat, momentum, and turbulence to an earlier-established depth of mixing. 


This happens when downward heat flux (buoyant damping) and dissipation effects 
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exceed the wind-stirring shear-production of turbulence, creating a stratification of the 
upper ocean. 

Thermal energy and mechanical energy received from the atmosphere not only 
control the local dynanucs, but the laver itself modulates the flux of this energy to the 
deeper water masses. Entrainment also converts some of the mean flow energy into 
turbulent energy, over and above the wind-stress production. 

Finally, substantial barotropic and baroclinic features, such as tidal motion 
and internal waves, can be linearly superimposed. The mean fields of concern are 
therefore the horizontally homogeneous components of the total fields. 

3. The model and its features not previously demonstrated 
he vertical and horizontal components of the turbulent kinetic energy (TKE) 
are determined implicitely, along with layer depth. mean momentum, and mean 
buovancy. Layer growth and retreat are predicted. 

Specific features of the TKE budget include mean turbulent field modeling of 
the dissipation term, the energy redistribution term, and the term for the convergence 
of buovancy flux at the stable interface as shown in a following section. Then an 
entrainment hvpothesis dependent upon the relative distribution of the TKE between 
horizontal and vertical components permits the closure of the system of equations. 

Wie amodel differs from earher models in the following ways. First. the 
amount of wind generated TKE to be used in mixing ts a function of the ratio of the 
mixed layer depth to the Obukhov mixing length L. Second, viscous dissipation is 
dependent on a local Rossby number. Finally, separate vertical and horizontal 
equations for TKE are used. permitting a more consistent interpretation of both 


entraining and retreating mixed layers. 


B. FORMULATION OF THE EQUATIONS USED IN THIS OBL MODEL 
1. Generalities 

The underlying principles employed in studving the mixed laver are the 
combined conservation of mass, momentum, thermal energy, and mechanical energy. 
Most of the physics behind such a one-dimensiona! model is based on the flux form of 
the Navier-Stokes equations with the Boussinesq approximation. The horizontal 
homogeneity mentioned previously permits the neglect of the horizontal gradients of 
the mean fields. 

Conservation of buovancy is employed as a generalization of the conservation 
of heat alone. The buovancy equation is generated from the heat and salt equations 


together with an equation of state, 


Cod 
oe) 


Pp = Pp, [l - a(8-8,) + B(s-sp)] . (3.1) 
and the definition for buoyancy, 

b = S(P5-P) MDG - (3.2) 
Where Py is a representative density at the time and location of consideration, g is 


acceleration due to gravity, @ and {} are the expansion coefficients for heat and salt. 


All variables are separaied into mean and fluctuating components : 


¢ temperature T=] T 

e salinity Se SO) oe 

2 presstine p=Ptp 

* velocities U = Cee Ua 


v= V-+Viv 
eWay. 
® buoyancy b= B+ & 
Subscript g denotes gecstrophic components. 
2. Nfean buoyancy and momentum equations 
The first law of thermodynamics for an incompressible fluid and the 
conservation of salt mass, neglecting molecular fluxes, lead to the mean buovancy 


equation : 
OB, dt = -db’w'/éz + agQ/ PoC, B89) 
By invoking the Boussinesq approximation, dropping the negligible viscous terms, 


assuming incompressibility, and subtracting the geostrophic equations from the total 


momentum equations, we obtain the mean momentum equations : 


OU, ot = [V - du’wéz (3.4) 
OV.dt = -fU - dv'w'/éz . (333) 


Integration of (3.35), (3.4) and (3.5) over the entrainment zone from z=-h-6 to z=-h 


leads to the jump conditions for the turbulent fluxes at the bottom of the mixed layer : 








-b'w(-h) = AB éhidt (5o0) 

-uw(-h) = AU éh’ét (ey) 

-vw(-h) = AV Ghidt (3.8) 
where 

b’ = agT’ - Bgs’ (3.9) 
and 

AB = ugAT - BgAS . (5a1,0) 


The assumption of vertical homogeneity in the mixed layer permits the integration of 
(3.3), (3.4) and (3.5) from z=-h-6 to z=0, including the effects of entrainment stresses 
(3.7) and (3.8), and entrainment buoyancy flux (3.6). This vields the bulk relationships 


for mean buovancy : 


hé<B>,ét + AAB chet = agQo/ PoC, - b'w'(0) (i) 
hé<U>/ét + AAU 6hiét = -fh<V> - u’w(0) (3.12) 
hé<V>-ét + AAV Ghét = -fh<U> - v'w(0) (3.13) 


Lore) 
Car 


Where 
Q, = Q. -(Q. + Q, +Q,) in Wm?. (3.14) 
0 S c h b 
The value of Q, is the net solar irradiance at the surface, minus the long wave back 
radiation, minus the sensible heat flux, minus the latent heat flux. 
The function /\ is the Heaviside step function : 
A= 1 if éh¥ct> 0, 


A=0 if chien = or 


The brackets < = denotes a vertical mean thi ouch tne micamay cm 


a) GZ (215) 
and denotes horizontal mean, 
Dee 
=— dv 16 
(F= SLL) dx ds (3.16) 


As the time step used in the model is one hour, the surface boundary conditions are 


prescrived hourly : 


“uw (0) = T,(t). Pp ( Slee) 

-yw(0) = T(t), Do {3.13 

“ow (0) = g{ Bs'w(0,t) - awT(0,t) | 3.19) 
With 

© = p.C,U,,° (3.20) 


Oe) 
ON 


where is the air densitv, C, a drag coefficient and U,,, the wind speed at 10 meters, 
4 ’ d S 10 P 


and with 
sw(0) = S(E - P) (3.21) 


where E-P = evaporation minus precipitation in cm/sec. 
3. Mfean turbulent kinetic energy equation 
Subtracting the scalar product of (u,v.w) with the mean momentum equations 


from the mean equation for the total mechanical energy vields the mean TKE 























equation: 
| GE =i 7 ee 
= uw— eee ve — + bw -—w(—+-—)]-e ~ 0 Spat 
20 lu 62 ee : at a yi oe) 
where (1) (IT) COD SN 
E=uet+ yet w?. (3.23) 
The budgets for the individual components of TKE can also be formed 
ae ae — au 6 fw “2 ‘Guo — —- 
= iy eee zs as See qe Oe 25 ues (3.24) 
Eat Cle Ol a PyPx 3 3 
l oe —— et Wy pov’ = 
Si —_--_ + ——-—- .@) L v 3 Zs 
2 at “Silay, ar ale eo | See oe oe) 
léw- — w wp oo pow ¢ ee - 
wae = bw (oe —) + —-—-= + Qhu’w (3526) 
ze Ot Po Po?Z = 
(V) (IV) 


where 22 = (Q2, ©,,©,) is the rotation vector for the earth. 
The time rate of change of TKE ts usually smaller than the other terms and 
may be neglected. The term (I) represents the rate of mechanical production and is the 


dominant source of TKE. It is the rate of conversion of mean to TKE by the down- 


oe 


gradicnt turbulent flux of momentum. The term (II) represents the buoyancy flux and 
can be either a source or a sink. This term can become an important source as in the 
case of strong convective cooling in the autumin. The term (II]) is the divergence of 
the turbulent flux of TKE or turbulent diffusion of itself. Locally, at the bottom of the 
laver during occasions of entrainment, a net convergence of flux of energy 1s necessary 
to maintain the downward buovancy flux for deepening the mixed laver. The term (1V) 
represents viscous dissipation. Because local isotropv is assumed for the dissipation 
range, ¢ 1s divided equally among the component budgets. The terms (V) are very 
important terms which sum to zero by continuity in the TKE cquation. They cause a 
redistribution of energy among uw, v* and w’*. The terms (VI) are also redistribution 
terms, but they are due to rotation of the carth. In this study, we will neglect them 
because of the usually short integral time scale in comparison with the time periods of 
INnterese 

Assumption of vertical homogencity in the mixed layer permits the vertical 
integration from z = -h-6 to z = 0 of (3.24)+(3.25) and (3.26). By this step we 


introduce a new variable h. the depth of the mixed layer. 


Up to that point. if we set q”, the horizontal component of TKE, such as 


C= (S27. 
we have SIX vatiablesen, <= Be = ie <we> and <q?> for five equations 
(3.11}, (3.12), G.13), vertical integral of (3.24)+ (3.25) and (3.26). Thereforeeamsaere 
equation 1s needed to close the svstem. Besides, we aiso need a suitable modeling of 
the different terms of the integrated TKE equation. Following Garwood’s arguments 
{Refs. 1,5]. a synthesis of these derivations 1s given in the next sections. 
4. Modeling of the different terms of the integrated TKE equation 
a. Shear production and turbulent diffusion 
The vertical integral of terms {1) and (II!) mav be combined to give the net 


“wind-generated” rate of production : 


Ge he (1) + (L11)] dz (3.28) 


Or 


Crd 
CO 


G = mu? + 12 ((AU)? + (AV)] dh, dt (3.29) 
where u™ ts the friction velocity defined as: 
" Sy gee CE — 
i Oy wy (08)? . C0) 


b. iNet buoyant damping 
The integral of (I]) over the mixed layer gives the net buoyant damping for 


mre Whole layer : 


B = fy (tl) dz Can) 
or 

B= SAB -b’w'(0)} + mC Se FQar | dz (3.32) 
which can be rewritten as 

B = 12h b’w‘(-h) - 1/2 h u*b* (3.33) 
wiiiere u*b* is the downward surface flux of buovancy : 

out = -b'w’(0) + Gg PoC, RQ{I + a1 = 27h) - apeph| . (3554) 
The radiation absorption Q(z) has been modeled as 

Q= YRQ ev? (oe) 


where ¥ 1s the extinction coefficient for net solar radiation, and RQ. is the short wave 
fraction of the net solar irradiance. This model assumes that the net long-wave solar 


radiation is absorbed at the surface, which leads to: 


-b’w'(0) = ag pyC, [(1-R)Q, - (Q.+Q,+Q,)] + BgS(P-E) (3.36) 


of 


where (1-R)Q. corresponds to the short-wave incoming radiation. On the other hand, 

the short-wave radiation penetrates below the surface where it is absorbed following 
the exponential decay (3.35). 
c. Viscous dissipation 

For a fully turbulent mixed layer, viscous dissipation of the turbulence 

occurs primarily in the small eddies which are locally isotropic. The merepameen 


dissipation is modelled as follows: 


_ an eee = “ae 
D = fg & &2 ees ee een (n= (Bs (3.37) 
or 
=e 322 a a = bee a 
D=m,<E=s 4) hme (3935) 
Where Rg = u*/hf is a Rossby number for the turbulent boundary layer. The first 


term in D arises from the fact that the time scale of the largest eddies meaner 


Ve 


proportional to the mixed layer depth divided bv the rms turbulent velocity <E> 2. 
The second one comes from the assumption that, in deeper boundary layers, planetary 
rotation {time scale 1/f) turns the mean shear direction with depth and thus influences 
the geometrical aspects of the integral scale. 

ad. Redistribution of TKE 


The vertical integral of the pressure redistribution term 





0: pou 
R. = —— q 3.39 
1 h-6 po@x, = ae 


i$ an important source or sink term for the individual TKE budgets, even though 
R,+R,+R,=0. The bulk formulation is 


R, = na 2 (216 ae si avce se (3.40) 


The rotational redistribution terms are assumed to be of higher order, and are 


neglected in this study. 
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Closure hypothesis 


Garwood [Ref. 1] achieves closure of the problem by formulating the following 


entrainment equation : 


« ~ Age Le - 
ch'ct = ee > US avalon 


e. 


SUMMARY OF MODELED EQUATIONS 


a ; : “ = Sy : 2 
Waemavecrsix Vanidoles th = B=, <U2)<V>, <w->, <q-> and a final set of 


Sie cquations : 


* mean momentum equations : 


hé<U>idt + AAU Cchiédt = 


ne Vet + AAV chet = -fh<U> - vwi0) 


* mean buoyancy equation : 


ime Be ot + AAB dhjét = 


-h<V> -u'w(0) 


U2Qq / P9C, - b’w'(0) 


—— 


® horizontal integrated TKE component : 


3 | (AU)? + (AV)* Gh 


h<q°>) = m,u*° ar 


. ay(< E> See ) ee a 2m 3 (<E—E> 
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@ vertical integrated TKE component : 


a 


ne NE Sea 1 h b’w'(-h) Slt Senin eae: m,(< E> - ee eo re 


ct 


, + 
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(3.44) 


‘mM, inv iS 


We 


(3.46) 


* entrainment equation : 


éhiét = m,<w?> @<E> / hAB. (3.47) 


The different constants used in the model are the following : 


Representative density of air 
Representative density of sea-water 
Coefficients of thermal expansion 
Gravitational acceleration 

Specific heat at constant pressure 
Drag coefficient 


Coriolis parameter depending on latitude 


We can tune the model by specifying the different parameters m,, m,, my. m, as well 


as the extinction coefficient for net solar radiation y, and the short-wave fraction of the 


solar fadiatione 


Initial profiles of temperature T(z) (°C vs cm) and salinity S(z) (g/kg vs cm) have 


to be provided. The depth increment in the model is 100 cm. Finally, as the time step 


is one hour, the following boundary conditions have to be prescribed hourly : 


- 


P 


Wind stress in dynes‘cm* 
Incoming solar radiation in Wim 
Back radiation in W/m" 

Latent neat tux 

Sensible heat flux 

Evaporation in mm/hour 


Precipitation in mm/hour 


The values of u* and u*b* can be computed using equations (3.17), (3/1S)ygmae 
(3:37) andi( 3.26), 


Py; COULPEING THE TWOMODELS 


mo )6UGGENERALITIES 

The goal of this research is to couple the Oceanic Boundary Layer (OBL) model 
derived by Garwood in 1977 with an acoustic ray tracing program allowing us to 
Miameessome eilects of the atmospheric factors on the oceanic environment and, 
therefore, acoustic propagation 1 the ocean. 

The two previous chapters gave us a theoretical approach concerning these two 
models. Now, we are going to highlight the coupling of these two models, the inputs 
that the coupled model needs, and the outputs we can get from it. The block diagram 
of Figure 4.1 depicts the over plan, as followed in the next sections. 

Initially, we enter the OBL model with a set of boundary conditions, with some 
initial conditions, and with a set of parameters including the time-step of integration of 
the OBL mode! and the time interval between each resulting ray trace. [rom the 
predicted output of temperature and salinity profiles, we compute a sound-velocity 
profile which allows us to trace the paths of acoustic rays according to a set of 
geometrical parameters such as the source depth, the maximum range, and the initial 


+ 
angle of a rav.~ 


B. INITIAL CONDITIONS 

We have to initialize the coupled model with a temperature profile, that 1s, 
temperature T in Celsius versus depth z with a one-meter increment and a salinity 
Peome, that is, salinity Sem oike wersus depth Z with a one-meter increment. These 
profiles can be obtained from elimatological data [Ref. 7], as the ones we used to 
analyze some applications of this coupled model, or from an XBT (expendable 
bathythermograph) and a salinity profile from climatology. We can also assume a 
constant salinity profile based on an average salinity in the area studied because of the 


small effect of salinity variations on sound-velocity computauon. 


“In all of the previous chapters, we wanted to be consistent with the coordinate 
systems used in each reference, that is, [Ref. 4] for Chapter II and [Refs. 1.5} for 
maiaprer iii. Im Chapter 11, y is the downward vertical axis and, in Chapter II1. z is 
the upward vertical axis. Since the coupling is done through the output T and S 
profiles of the OBL model feeding the input profiles of the acoustic model. this 
notation conflict does not affect the coupling by itself, nor the codes of the coupied 
programnis. 
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Figure 4.1 Coupling the OBL model with the ray tracing program. 
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C. BOUNDARY CONDITIONS 

The boundary conditions or atmospheric forcing factors we have to provide are 
listed below. 

The hourly solar radiation flux Q. has to be given in W/m versus time in hours. 


The climatology gives us generally a daily average Q In our study we will simulate 


savg 
a diurnal cycle for Q. based on the following formula : 


Pee, sim l2) if O<t< 12 hours, (4.1) 


ee 0 ee Ours. 


as depicted in Figure 4.2. 





Figure 4.2 Simulated diurnal cycle. 


A straight-forward integration leads to the following relation : 


Oe eax ee Seen : (4.2) 


Simulation of clouds is done by allowing only a part of this Q. as input solar radiation 


(20 to 50% of Q. for example) at the sea surface. 
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The values of the latent flux of heat Q., the sensible flux of heat Q, and the back 
radiation Q, have to be provided on an hourly basis. In the simulations of the next 
chapter, we will pick some average values of the sum Q.+Q,+Q, from the 
climatology, depending on what time of the year We are working with. 

The evaporation rate EV in mm/hour is also evaluated from climatological data, 
and precipitation rate PR in mm/hour can be simulated. Heavy rain cases will be 


considered also because of some interesting effects. 





TOAD ie 2 
WIND SPEED AND WIND STRESS CORRESPONDie ies 


Wind Speed Wind Stress 

knots nis g/cm/sec” 

5 rae 0.106 

10 5 0.423 

es Tes USS 

20 10 1.694 
| 25 Zo 2.647 

30 15 Sued! 

35 Li 5.188 
40 20 6.776 


2 


The wind stress t in g’em/‘sec” has to be provided hourly also. “igen 


simulations, we will consider both constant strong wind as well as light wind 
conditions. Based on equation (4.3) 


(4.3) 


5 5 
Where an air density p, of 1.21 kg/m” and a drag coefficient C, of 1.4x10° have been 


chosen. The correspondence between wind speed and wind stress is given in Table 2. 
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D. OTHER INPUT PARAMETERS 

The other parameters we have to input in the coupled model are the latitude of 
the area studied, the number of davs NDAYS during which we want to integrate the 
model, the interval of integration used by the model (IDIFF = lhour in this study), and 
the frequency TI in hours with which we want to output some graphic results from the 


acoustic ray tracing program. 


Pee OULPUTS OF THE OBU MODEL 

Poi tic Ob leinodel We-cam eeu tie nourly values of the predicted depth, 
temperature and salinity of the mixed layer, and we are able to follow the evolution of 
these different parameters on tabulated outputs. 

The link with the acoustic model is done through the hourly computed profiles 
T{z) and S(z). Depending on the frequency TI, these profiles (renamed T(v) and S(v) 
for notational consistency as mentioned in a previous note) are used to compute the 


sound-speed profile c(v) by a special subroutine SVEL. 


Fe. SOUND-SPEED COMPUTATION 

In the subroutine computing the sound-speed, we used the equation derived by 
Chen and Millero {Ref. 8]. The sound-speed obtained agrees with the data of Del 
Grosso at one atmosphere [Ref. 9] and with the data of Wilson at high pressures 
[Ref. 10}. Other formulas for sound-speed are mentioned by Urick [Ref. 11]. but the 
One we used is miost often suggested for oceanographic calculations. Further more, 
comparisons of graphic outputs do not reveal any significant differences between the 


results obtained by using these various equations. 


G. RAY TRACING 

The computed sound-speed profile c(y) 1s used as a main input in the acoustic 
model. Depending on the case we want to study, a set of geometrical parameters has 
to be provided such as the source depth YO, the maximum range RANGE, and the 
maxinium depth N. The bundle of rays to be plotted is defined bv the number of rays 
M, the initial angle BO of the upper ray of the bundle, and the angle increment DELB 
between two rays. Different equations for defining a bundle of ravs can be set inside 
the subroutine. Thus, it is always possible to pick up a specific bundle of rays to 
highlight some effects of the atmospheric forcing on some precise type of rays and 


Gemonstrate some acoustic energy redistribution. 
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Since, in this research, we want to give an example of the application of the use 
of this coupled model, we will choose generally a ray-bundle width of 5 to 8 degrees 
which should correspond to an active sonar operating at a frequency in the range of 5 
Rize 

Finally, the main output of the coupled model is a plot of the sound-speed profile 


aside an acoustic ray tracing for a given set of parameters defined previously. 


H. NOTE ON THE ALGORITHM USED IN THE ACOUSTIC MODEL 
Usually, in most acoustic models, the input sound-speed profile is approximated 


by a piecewise linear curve as depicted in Figure 4.3. 
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Figure 4.3 Piecewise linear approximation 
of a typical sound-speed profile. 


Thus, the profile can be modeled mathematically by the following equation: 


CY) = Cae 2 vec , Yep SY SY; (4.4) 


and is composed of several distinct straight line segments, each with its own constant 
gradient g.. In each such delimited slab i, because of the constant gradient g,. the ray 
path is an arc of a circle for which an equation is easy to compute and fast to plot. 
Since a profile is usually divided into a few discrete slabs (for example, 10), a more 
computationally efficient subroutine could be derived from the one discussed in 
Chapter II. This subroutine would invoke gradient sign checking and initial angle 
comparison to the horizontal only when a ray crosses a boundary between two slabs. 

In our study, this checking and comparison are done at each depth increment (1 
m) as well as the plotting of a new point belonging to the ray path. It is certainly a 
much slower algorithm than the previous one, but we had to consider that, even if the 
output profiles T(y) and S(v) of the Garwood OBL model can be approximated as 
piecewise linear, the output sound-speed profile c(y) cannot be considered as piecewise 
linear. This is because c(y) is a highly non-linear function of T. S, and y. This 
characteristic can be demonstrated by using a simple equation like the one derived by 
Moppens (Ref. 12] 


c(T,S,y) = 1449.05 + 4.57T - 0.0521T* + 0.00023T? (4.5) 
+ (1.333 - O.0L26T)(S - 35) + 0.0163y. 


menis may be rewritten as 

c=a+t bl + ct? + dT? + eS + FITS + gy. (4.6) 
Differentiation of c versus depth v vields 

de dy = g + (b+ 2cT+3dT7+fS) dT/dy + (e+fS) dS'dy. (4.7) 
Hence, in a particular slab i, even if (dT, dy), and (dS’dy), are constant, ipacall berceen 
from (4.7) that (de/dy), is not a constant because T and S are not constant in the slab, 
but are linear functions of depth y. 

Finally, the acoustic model developed in Chapter II permits the use of ABT 


profiles discretized with a depth increment of one meter or more. Use of climatology is 


often unrealistic due to multiple averaging and too coarse vertical resolution. 
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V. APPLICATION OF THE COUPLED MODEL TO A SPECIFIC AREA 


A. GENERALITIES 

In order to investigate the typical results we can obtain from this coupled-model, 
we shall simulate some cases in the Western Mediterranean Sea. The chart presented 
in Figure 5.1 shows the site to be studied: @=42.5°N, A =007.5°E. 
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Figure 5.1 The northern region of the Algero-Provencal Basin. 


The Mediterranean Sea is one of the major, deep closed seas of the world. The 
water mass interacts with the open ocean only through the Strait of Gibraltar, which 
has a sill depth (300 m) much shallower than the average depth of the rest of these 
closed waters. The water below this sill remains almost homogeneous in temperature 


and salinity. The western Mediterranean is an area having very specific oceanographic 
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and atmospheric features. Figure 5.1 depicts the north part of the Algero-Provencal 
Basin, which is one of the two large deep basins of the Western Mediterranean. The 
basin bottom contains wide, flat abyssal plains, averaging 2700 meters in depth. The 
Water mass 1s Clearly bounded by land and islands, and this sea is well ventilated due to 
the quantity of deep water formed at the surface by winter cooling. As described in 
(Ref. 13], during winter, the Alpine relief brings intrusions of polar continental air 
masses into contact with the surface waters of the northern part of the Mediterranean. 
Strong, cold, dry, continental winds (mistral) can blow for a few days, cooling the 
surface water and leading to a vigorous convective mixing. Such a strong sinking 
motion of surface water can mix and cool uniformiuy as much as the upper thousand 
meters of the water column. The sea-state is often characterized by a short wavelength 
swell which causes a very rough sea during storms. Another interesting result of the 
interaction between the atmosphere and the sea is that the Mediterranean water 1s 
characterized by high temperatures (13°C at 1000m) and high salinities (38.4-38.45 
g/kg) which are surpassed only in the Red Sea. This 1s due to the large heat input and 


an excess of evaporation over fresh water input (precipitation and river contribution). 
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Figure 5.2 Seasonal variation of the temperature profile 
in the northern part of the Western Mediterranean. 


Figure 5.2 from (Ref. 14] gives an example of the seasonal variations observed in the 
Ligurian Sea, that is, the northern part of the Algero-Provencal Basin. It can be 
observed that, in this region, the water column is completely isothermal in winter. A 


warmer surface layer evolves progressively during springtime. In autumn, it decreases 
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in temperature but increases in depth until winter conditions are reached again. The 
summer temperature profile 1s characterized by a top layer about 20 to 40 meters thick 
overlving a thermocline with large temperature gradients (up to | to 2°C/m ). 

The Mediterranean Sea is also of interest for the numerous naval operations 
conducted in that area. For this reason, it is of interest to studv some of the features 
of the acoustic propagation that are related to atmospheric factors. For example, in 
this region, the heating by the sun of the upper layers of water, together with an 
absence of mixing by the wind, causes a strong near-surface negative scund-speed 
gradient to develop during the spring and summer months. This thermocline overlies 
isothermal water at greater depths. The result 1s a strong shallow internal sound 
channel (SOFAR channel: sound fixing and ranging) with its axial depth near 150 
meters. During the autumn, the profile returns to its winter time conditions, with 
isothermal water and positive sound-speed gradient extending to the surface of the 
water columin. leading to a half-duct tvpe of propagation. During the summer season 
the near-surface negative gradient and the resulting strong downward refraction greatlv 
limit the detection range of surface ship hull-mounted sonars. By way of 
compensation, the summer time channel in the Mediterranean produces convergence 
zones for a near-surface source similar to the deep-ocean sound channel, although the 
range of the zone is much less because of the smaller vertical extent Of (neve 
compared to the open ocean. Therefore, from the point of view of sound propagation, 
the important factor here is the existence of a predominantly isothermal and isohaline 
water mass at depth, which results in a sound-speed profile characterized by a deep 
constant-gradient section. Only the upper section of the profile will contain 


irregularities as depicted in Figures 5.3 and 5.4. See [Refs. 11,14] for more details. 


B. SUMMARY OF THE DIFFERENT PARAMETERS USED IN THIS STUDY 
ln the simulations discussed in the next sections, we used the parameters listed in 
Table 3. Four different cases are to be analyzed, corresponding to the months of 
December, February. June. and September. The boundary conditions are assumed to 
be climatological and were obtained from [Refs. 16,17,18.19]. A precipitation rate of 2 
mm hour over 12 hours will be assumed for heavy rain cases. Values of 20 to 50% of 
Q. will be used to simulate the reduction of radiation incident on the sea surface due to 
the presence of clouds. The initial profiles of temperature and salinity come from the 
climatology used by the Fleet Numerical Oceanographic Center (F NOC). Finally, we 


will avoid the complication of bottom reflection by assuming an infinitely decp ocean. 
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Figure 5.3. SSP for the Ligurian Sea, Western Mediterranean 
February to July from left to right. 
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Figure 5.4 SSP for the Ligurian Sea, Western Mediterranean 
August to January from right to left. 
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1540.0 


& DEGEYPRER CASE 
1. Source at 10 meters 
This case simulates a surface-ship hull-mounted sonar and we will consider 


strong wind (30 Knots) as well as light wind (5 Knots) conditions. 


| | 
| | 


| er ier 
Pier eer I IATOLOGY FOR THE MEDITERRANEAN 
(WM?) 
December February June September 
Q. 130 190 250 190 
One 400 600 780 600 
Or QO. 190 120 80 120 
Q,+Q,7Q, 280 190 130 190 
OF. -150 v 120 
EV(mm/hr) 0.3 0.2 0.1 0.2 





a. Strong wind conditions (30 knots) 


Iie evolution Of @itesmmxed layer depen h for different cases is shown in 
Table 4. In the case of no clouds and no rain, we find that the mixed layer depth 
increases significantly due to the combined effects of the strong winter ume surface 
Boone (Q_..= -150W'm") plus the effect of the wind stirring. 

Figure 5.5 shows the effect of this mixing and cooling on the acoustic ray 
paths. The first graph shows the acoustic propagation at time = OQ hour for a 
beamwidth @ = 6°, and the second graph shows the propagation after 24 hours of 
strong wind. By counting the rays, we find that 39% of the acoustic energy has been 
redistributed. More ravs are filling a deeper mixed layer, and the convergence zone 
1G7) has weakened. There is no effect on the CZ range which remains between 18 
and 20 km. 

Figure 5.6 shows, with more detail, the transformation of a bundle of 
refracted-surface reflected (R-SR) ravs (initial angles between 87° and 89.3") after 24 


hours of strong wind. For a total of 24 R-SR rays, 13 rays are now trapped in the 


>The term convergence zone (CZ) is applied to a phenomenon of focalization of 
the underwater sound in annular domains. 
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Figure 5.5 December, strong wind, source at 10 meters: t 
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sformation of R-SR ray 
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Figure 5.6 December, tran 


mixed layer, increasing the possibility of detection on a shallow target for a surface- 


ship sonar. On the second graph of Figure 5.6, the weakening of the CZ is quite 


apparent. 
TABLE 4 
EVOLUWMN OF THE MINEDSEA aR air 
IN DECEMBER WITH A STRONG WIND 
no rain no rain no rain heavy rain for 12 hrs 
| no clouds 50% Q. 20% Q. 20% Q. 
t(hr) h(m) h(m) h(m) h(m) 
0 37.0 37.0 BO 37.0 
| he ol SZ Sol 50.6 
| 24 60.3 Gile2 Oley 59.4 
| 36 65.7 67.0 Gil 65.7 
| 48 70.8 72.0 m8 70.9 
| 60 ae 76.0 76,9 ol 
Tt 782 WS 80.4 78.7 


From Table 4, we can see that adding clouds to the simulation increases 
the surface cooling and thus accelerates a somewhat the deepening of the mixed laver 
due to the added production of turbulence by wind simmme. The eliecmman 
superimposing heavy rain for 12 hours is a slight damping of the mixed layer deepening 
rate. The main conclusion of these last test cases is that, with strong winds, the overall 
effect of clouds and rain is not too significant and the general shape of the acoustic ray 
trace does nel change: 

b. Light wind conditions (3 knots) 

Light wind conditions allow us to analvze the interesting effect of heavy 
rain on the acoustic propagation. The evolution of the mixed layer depth is shown in 
Tables. 

Figures 5.7 and 5.8 show the effect of a heavy rain on a winter time sound- 
speed profile. 


The gradient of the sound-speed versus depth can be expressed as : 
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Figure 5.7 December, light wind, heavy rain for 
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Figure 5.8 December, light wind, heavy rain for 1|2hrs, 


dejdv = dciéT dT/dv +0c;dS dSidv + écidp dp/dy (5.1) 


where the average values of the partial derivatives are the following : 


oc.6T = 4.0 misec’°C , (ene) 

ecloS = 1.2 misec'g/ks (5.3) 
and 

ociép = 0.016 mfsec/m. (5.4) 


TUNG es | 


Ey OU TON Opel NED WAN ERSDEPTH IN DECEMBER | 
Wea GT ese SI HEAVY RAIN FOR 12 HRS | 


t(hr) h(m)} t(hr) h(m) 
0 0) 24 39.9 
6 deg 36 43.0 
12 10.4 4§ 46.4 
hs 14.4 60 49.3 
14 34.5 72 Se 





The temperature term of the gradient dc/dy is dominant, and the salinity 
term is usually so small that it is very often neglected. In the case of heavy rain, 
however, the term 6c/cS dS,dy is not negligible and can be very important and even 
dominant. In the mixed layer, dT/dvy=0. Thus if heavy rain occurs, we get 
dS,;dy>>0, and dce'dy can be positive and much greater than usual. This yields a 


sharp “inversion” as depicted by the sound-speed profile of Figures 5.7 and 5.8. 
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Another way to examine the effect of the rain is to consider the surface 
buoyancy tlux By : 


hoe 
— 


By = agQ,../ PC, - BeS(EV - PR). (5.5) 


In winter time, the first term of B) is on average strongly negative. This corresponds 
to an upward buovancy flux. During heavv rain fall, EV-PR is negative, leading to a 
positive second term in By. This fact can be verified by investigating the last two 
columins of Table 4. Rain fall damps the upward buoyancy flux and, thus, decreases the 
nuxed laver deepening. 

In Figures 5.7 and 5.8, we see the effect of heavy rain fall on the acoustic 
propagation at time 0, 6, 12, and 14 hours. In these graphs, the beamwidth has been 
set to 5°. After 6 hours of heavy rain, a high density of R-SR rays have been trapped 
in the mixed layer. After 12 hours, the effect is even much more important. Almost all 
the rays are trapped in a shallow mixed layer, and the convergence zone has 
disappeared. At that time, the only chance for a surface-ship sonar to get anv 
detection 1s for a shallow target. 

However, this effect does not last very long after the rain stops. We can 
see in Table 5 that, in onlv two hours, the mixed laver depth almost returns to its 
Initial value. We can also notice in Figures 5.7 and 5.8 that, only two hours after the 
end of the rain, the acoustic propagation is almost the same as for the initial case at 
t=O. This 1s due to the strong Winter time and night time surface cooling, even though 
the wind is light. This fact also demonstrates the strong effect of surface cooling on 
the muxed laver deepening in winter time. 

2. Source at 40 meters 
Now, We are going to studv the case where the source is just below the initial 
MuNedslamer (a tan}: 
a. Strong wind conditions (30 knots) 

In the case of no clouds and no rain, Figure 5.9 shows the evolution of the 
acoustic propagation over one day of strong wind for a beamwidth 0 of 6°. 

Initially, at time t=O, all the ravs are R-SR, and we have a perfect case of 
convergence zone (CZ) detection with a strong focusing effect at the smaller range of 
the CZ (around 15km). After 24 hours of strong wind, the deepening effect of the 


nuxed layer redistributes 42% of the acoustic energy into a mixed layer type 
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Figure 5.9 December, strong wind, source at 40 meters: t= 0,24 hrs, @=6°. 
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propagation. The CZ range stavs the same, but the CZ is weakened and the focusing 
effects ereatvereduced: 

Figure 5.9 shows that the detection ability of a towed sonar at 40 m is 
much improved because of the deepening of the mixed layer mainly due to the wind 
stirring. For the purpose of highlighting this effect, Figure 5.10 extracts only the 
influenced ray bundle, that is, the bundle with a beamwidth of 2.5°. All the ravs of this 
bundle are transformed from RR (refracted-refracted) type to mixed laver trapped. 

If we would have added clouds to the simulation, with or without rain, we 
would have found the same general shape for the acoustic ray trace. The fact confirms 
that, in winter time, the wind mixing has a dominating effect on the acoustic 
propagation. 

b. Light wind conditions (5 knots) 

For a source at 40 meters, even with light wind conditions, acoustic 
propagation 1s not very sensitive to heavy rainfall. Figure 5.11 shows the ray trace at 
time 0 and after 12 hours of rain. The difference in shape between the two graphs is 
not very great. Thus, the effect of heavy rain is to be considered significant, even with 
light wind conditions (which is not usually the case in winter time during storms 
passing by), for only a shallow transmitter, as seen in the previous section. 

3. Source at 100 meters 

If the source is sufficiently deep relative to the mixed layer Gepreine 
atinospheric forcing has no effect on the acoustic propagation. Adding clouds and rain 
will not change the general shape of the acoustic fray paths. Pieure 5.1/2 depicieaine 
case of strong wind conditions with no clouds and no rain at t=0 and 24 hours. Fora 
beamwidth of 6°, the RR rays stay trapped in the shallow SOFAR channel and do not 
enter the mixed laver. This is because of the strong thermocline (strong negative 
gradient dc,dy) underlying the mixed laver. For rays to enter the mixed layer, we 
would need steep initial angles at the source, which is not realistic for the usual sonar 
beaniwidth. 

In Figure 5.13, with a source at 80 meters, an interesting effect on (iewn 
paths is depicted. The rays almost appear to be reflected on the the sharp 
discontinuity created at the top of the thermocline because of the strong wind mixing. 

4. Source on the SOFAR axis (150 meters) 
Figure 5.14 shows the acoustic ray paths obtained when the depth of the 


source is set to 150 meters, that is, the optimum depth for the internal sound-channel 
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Figure 5.10 December, strong wind, source at 40 meters: t= 0,24 hrs, 0= 
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Figure 5.12 December, strong wind, source at 100 meters: t=0,24 hrs, 8=6°. 
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Figure 5.13 December, strong wind, source at 80 meters: t=0,24 hrs, b= 91° to 93°, 
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Figure 5.14 December, strong wind, source at 150 m: t=0,48 hrs, 8= 9°. 


propagation. Strong wind conditions were simulated for 48 hours and the beamwidth 
was set to 9°. The interesting acoustic propagation is due to the sharp positive and 
negative gradients dc/dy under and overlying the depth of the SOFAR axis. As 
mentioned previously, the source 1s deep enough for the acoustic propagation not to be 
influenced by the atmospheric forcing. A weak CZ appears at ranges between 15 to 18 
km, and these ranges are influenced almost not at all by the boundary conditions. 

In that case, if we superimpose the second graphs of Figures 5.5 and 5.14, we 
can see the combination of the ray traces obtained with a source at 10 meters and a 
source at 150 meters. The water mass is almost completely filled with acoustic energy, 
leading to a high probabilitv of detection of a target whose depth would be between 0 
and 300 meters. 

5. Source at 500 meters 

Figure 5.15 displays the case where the wind has been blowing for 24 hours 
With no clouds and no rain. We would have obtained the same general shape by 
varying the wind and (or) adding clouds and rain. As in the two previous cases, the 
source is deep enough for the acoustic propagation not to be influenced by the 
atmospheric factors. In any case, these graphs were provided to show the excellent 
shape of the acoustic rav trace in the case of a perfect surface reflection. 

6. Conclusions for December 

On short time scale, the atmospheric conditions have an effect on the acoustic 
propagation onlv for a shallow source {10 to 40 meters in our simulations), that is, for 
a source in the initial mixed laver or just below. 

For such a shallow source, the wind has a dominating effect over rain and 
clouds. Mainly, this effect 1s the transformation of R-SR rays to mixed laver trapped 
rays and the weakening of the convergence zone. 

The effect of the surface cooling is also important in winter time. The net 
upward heat flux was set to -150 W,m? in our simulations. This leads to an average 
deepening of the mixed laver untill the month of February when the entire water 
column is mixed. This was depicted in the light wind conditions cases, and it vields a 
haif-duct type of acoustic propagation. 

The effect of heavy rain has been also depicted for light wind conditions and 
leads to a strong shallow inversion, trapping most of the ravs emanating from a 


shallow source in a shallow surface duct. 





0°004 
. 0°009 

= < = == = = =a : 0°00s 
‘ ve U, , 0°00 

) / 0°00z 

x «s a AS j 0°00! 
leltentonre 8 


0°O£ 0°Sz 0°02 o°sI 0°! 0°S 0°70 ardisi 0°20SI 
(WH) GONWY (S/W) dSS 





0002 
0° 009 
=e 000s 
oo 0°O0F 
RK 0° 00s 
RK ; 
SR ROX 0°00z 
OY : | 
G x Rey 0°00! 
SOX Kaen 
K KROQ KOS ee ae 0°0 
Oot 0°S2 0°02 0°St 0°01 0°s 00 OrZISI 0° £0S1 


(WH) GONBY (S/W) dSS 


6°. 


(W) HLd30 
g , source at 500 m: t=0,24 hrs, 8 
ol 


strong wind 


(W) HLd30 
December, 


Figure 5.15 


D. FEBRUARY CASE 
1. Source at 10 meters 
a. Strong wind conditions (30 knots) 

Table 6 shows the evolution of the mixed layer depth for different cloud 
and rain conditions. In all the cases, the mixed layer deepens mainly because of the 
wind stirring and, as we saw previously in the December case, adding clouds to tne 
simulation increases the deepening of the mixed layer, while the rain reduces slightly 


this effect. 





TABEEsc 


BVOLEUTION OF THE Ai Se DE ear ii 
IN FEBRUARY WITH A STRONG WIND 


no rain no rain no rain heavy rain for 12 hrs 
no clouds 50% Q. 20% Q. 20% Q. 
t(hr) h(m) h(m) h(m) h(m) 
0 70.0 70.0 70.0 70.0 
| i 77.6 TS 80.5 76.8 
ape Sons 87.3 88.3 84.9 
36 89.1 92.3 94.1 90.8 
48 94.5 97.7 99.5 96.3 
60 96.9 LOS 103.8 100.7 
72 2 105.6 108.0 104.7 





The initial mixed layer depth extracted bv the OBL model is 70 meters. In 


this model, the criterion on which is based the choice of this depth is the following : 
Ab) = |agAT - BgAS| > 0.04 cm’s? (5.6) 

where the symbol A represents the variation of buoyancy, temperature, and salinity for 

an increment of one meter depth. A higher value for the criterion (0.05 for instance) 


would have increased all the values of the mixed layer depth, but the evolution would 


have had the same trend as the case depicted in Table 6. 
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However, independent of the choice of the criterion for the mixed layer 
depth, the sound-speed profile of Figure 5.16 shows that the water column is mixed 
almost from the top to the bottom with an average sound-speed gradient de dv of 0.02 
sec). Thus, in February, the main characteristic of the acoustic propagation ts a half- 
duct type of propagation. In the case illustrated by Figure 5.16, the beamwidth has 
been set to 12°, and the two acoustic ray traces show the situation at time 0 and 24 
hours. The wind does not affect the acoustic propagation since the water mass is 
already almost completly mixed at t=Q. Nearly identical rav path shapes would have 
been obtained by assuming clouds and rain. 

b. Light wind conditions (5-10 knots) 

Table 7, for hght wind, no clouds, and no rain conditions, highlights the 
diurnal variabilitv of the mixed layer depth and illustrates how the surface waters of the 
sea Warm during the daytime hours on a sunny day, shallowing the mixed layer depth. 
At night, the deepening and cooling mixed layer returns to a state nearly the same as 
the initial one. These diurnal changes have a profound effect on sound transmission 
for a surface-ship sonar as depicted in Figures 5.17 and 5.18. 

[In Figure 5.17, we simulate a light wind of 5 knots with no clouds and no 
rain. The ray traces are shown at time instant 0 and 12 hours for a beamwidth of 8°. 
At the end of the daytime (t=12 hours), the general shape of the acoustic ray trace is 
almost the same as the one at the end of the might time (t=0 hour), except for the 
ranges less than 6 kin tn shatlow water. This is depicted tn Figure 5.18 at t=§& hours, 
Heat is, im the middle of the afternoon, for a beamwidth of 6°. This reduction of 
suriace-ship echo ranging on a snallow target is a phenomenon often called the 
“afternoon effect.” 

Next, heavy rain is included for 12 hours with hght wind conditions. The 
effect is depicted at time 0 and 12 hours in Figure 5.19, where the wind speed was set 
equal 10 knots and the beamwidth to 6°. At the initial time, the acoustic energy 
propagates following a half-duct. At time 12 hours, we find that 62% of the ravs are 
trapped in a shallow surface duct because of the strong discontinuity in the sound- 
speed profile due to the decrease of salinitv in the upper laver of the sea caused by a 
weak wind stirring and damping of the turbulence by the strong downward surface 
buovancy flux associated with the rainfall. However, the overall shape of the acoustic 
ray trace is not changed significantly. There will be an improved possibilitv of 
detection for a shallow target. If the wind continues to blow at 10 knots, the effect 


remains even 12 hours after the rain has stopped. 
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Figure 5.16 February, strong wind, source at 10 meters: t 
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Figure 5.18 February, light wind, source at 10 meters: t=0,8 hrs, 6 = 6°. 
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2. Source at 200 meters 

In this case, we simulate the performance of a sonar towed at a depth of 200 
meters. In order to demonstrate that the atmospheric forcing has no effect on sound 
transmission from a deep source, Figure 5.20 shows the case with a heavy rain for 12 
hours together with a light wind of 10 knots. We can see that rain does not influence 
the acoustic propagation. Whatever the atmospheric conditions in Februarv, a deep 
towed sonar 1s able to detect any target moving between 0 and 300 meters. Thus, the 
half-duct type of propagation is one of the best configurations for underwater detection 


by deep transmitters: 












Pepe 7 


EVOLU TIGN OF THE MIXER aie aioe ee rants 
IN FEBRUARY WITH A LIGHT WIND 





wind 5 knots wind 5 knots wind 10 Knots 
no rain heavy rain heavy rain 
for 12 hrs for 12 hrs 
no clouds 20% c 20% Q, 
t(hr) h(m) t(hr) h(m) t(hr) h(m) 


0 0) 








3. Conclusions for February 


As in December, the atmospheric conditions do not influence the acoustic 
propagation for rays emanating from a deep towed sonar. 
Compared to December, the SOFAR channel has disappeared, and there is a 


completely asymmetrical, surface-bounded channel (or half-duct) inside which the 
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Figure 5.20 February, light wind, heavy rain, source at 200 meters: t 
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propagation is made of R-SR rays only. Because of this half-duct type of propagation, 
February is one of the best months in the Mediterranean for underwater detection by 
both deep and shallow transmitters. This generalisation not take in consideration the 
sea state Which can generate noise and reflection loss when the sea 1s rough, because 
our acoustic model assumes a perfect surface reflection. 

With light wind conditions, surface-ship echo ranging can be poorer in the 
afternoon on a shallow target because of the “afternoon effect” due to surface water 


heating. 





Figure 5.21 Variation of sound-channel width d for similar temperature 
changes in (a) the Mediterranean and (b) the Atlantic. 


EF. JUNE CASE 

l-igures 5.3 and 5.4 show that the shape of the sound-speed profile stays 
approximatively the same during the entire sumimer season, that is, from June to 
September. Therefore, in this section, we will just simulate the atmospheric conditions 
and the acoustic propagation in June. All of our observations will be applicable to the 
other months of summer time even though the atmospheric factors change somewhat, 
as illustrated in Table 3. 

In summer, so much heat is added to the upper layer of the Mediterranean Sea 


that negative gradients in temperature and, thus, in sound-speed develop in the upper 


SO 


laver of the sea. This 1s underlain by isothermal water at a relatively shallow depth. 
The consequence is that a sound channel exists in summer in the Mediterranean and ts 
characterized as follows: 

The sound channel 1s, in all cases, strongly asymmetrical. The sound-speed 
profile contains a sharp discontinuity between sound-speed gradients of completely 


different orders of magnitude: 0.017 sec”! 


-| 


Om (hem@ ver side, due to the pressure effect, 


and -1 to -5 sec”* on the upper side in the thermocline. 


eB bes 


Ey Oeetereor fine VlXED LAYER DEPTH 
Pe ie fli STRONG WINDS 


| 

| 

no rain no rain no rain heavy rain for 12 hrs 

no clouds 50% Q. 20% Q. 20% Q. 
t(hr) h({m) h({m) h(m) h(m) 

| 0 2.0 2.0 2.0 2.0 

| 12 20.3 20.8 Ziel 20.3 

24 De 26.2 26.5 Dore 

: 36 29.4 30.3 30.8 30.0 

48 32.9 oo 34.5 33.6 
60 Ehie 36.9 Bd 36.8 | 





5 
| ie 38.3 39,8 40.7 39.8 
| 


The sound channel axis is verv close to the surface, and the curvature of the 
sound-speed profile near this minimum is very great. 

The width of the sound channel is smaller than in the Atlantic and is extremely 
Pamdple with the temperature of the surface water, as depicted in Figure 5.21. The 
usual width of the sound channel in summer is on the order of 1400 to 1800 meters. 
The depth of the Mediterranean is sufficient for this channel not to be intercepted bv 


the sea bottom in most cases. 


81 


[. Source at 10 meters 
a. Strong wind conditions (30 knots) 

The evolution of the mixed layer depth h for different cases is Siow 
Table 8. In the case of no clouds and no rain, we find that the mixed laver depth 
increases significantly during the first 12 hours because of the strong wind mixing in 
spite of the strong summer time heating. As seen in the previous sections, for short 
timie scales, the wind is the dominating factor influencing both the mixed layer depth 
and the acoustic propagation. 

Figure 5.22 shows the acoustic ray trace at time 0 for a beamwidth 0= 8°. 
After 24 hours of 30 knot winds, the mixed laver deepens to 25,7 mieters. Evenmnermc 
convergence zone weakens slightly, the range of the CZ is not affected. However, the 
acoustic propagation changes greatly for the less steep rays. Thirty percent of the rays 
are now trapped in a shallow mixed layer, improving slightly the probability of 
detection using a surface-ship sonar. However, the unfavorable conditions for 


underwater detection for a shallow source 1s well depicted in Figure 5.23. 


TABI ES 


MAXIMUM AND MINIMUM MU DEPTHS ASipetEMrPER he 

IN JUNE FOR LIGHT WIND CONDITIONS Wit @ Geo ee. 
| 
| 
| 


day 2 9a.m. h=4.8m T= 20.13°C 
Sijolio0 les 1 Uae T=21.3PC 
dav 3 9a.m. h=5.4m T=20:15°@ 
5p.m. h= 1.1m T=2i.33 € 
day 4 9a.m. h= 5.6m = 20ee tC 


| 
| Sperm. h=1.lm l=? 
| 
| 


b. Light wind conditions (5 knots) 
Table 9 shows the evolution of the maximum and minimum mixed laver 
depths and high temperatures rise during light wind conditions with no clouds. The 
mixed layer responds to the diurnal cycle of solar radiation, shallowing in daytime 


because of solar heating, and deepening at night because of surface cooling. We notice 
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Figure 5.22 June, strong wind, source at 10m: t=0,24 hrs, @= 8°. 
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Figure 5.23 June, strong wind, source at 10m: t=0,24 hrs, 8= 8°. 
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that the general trend of the minimum and maximum temperature is to increase, 
leading progressively to the extreme sound-speed profile of August (see Figure 5.4). 
Nethertheless, the mixed layer stays sufficiently shallow not to influence the acoustic 
propagation, even for a shallow transnutter. 

Adding clouds to the simulation (20% Q_) only deepens the mixed layer to 
9 meters after 48 hours, which has almost no effect on the acoustic propagation. A 
longer period of cloud covering would not be a realistic simulation for summertime in 
the Mediterranean. 

Finally, heavy rain keeps the nuxed layer shallow (around 1.5 meters), and 
sound propagation is not influenced. 

2. Source on the SOFAR axis (90 meters) 

Figure 5.24 simulates the case where the source is on the shallow sound- 
channel axis characterizing the Mediterranean summer time. The beamwidth is set to 
8°. The strong asymmetry of the sound channel is readily apparent. As previously 
mentioned, the wind does not influence the sound channel propagation for 8=8°. The 
fueseat the Source are not sufficiently steep to be affected and, even if the sound 
channel width decreases, the propagation does not change. 

However, in order to study the effect due to the decrease of the sound channel 
width, Figure 5.25 shows the ray traces for steep initial angles B= 81°, 82°, 83°, 97°. 
p . and 99 at time 0 and after 48 hours of strong winds (30 knots). The channel 
width decreases from 1400 to 900 meters. The steepest rays (f= 81°, 82°, 98°, and 99°) 
are transformed from RR type to R-SR type, but the influence of this phenomenon on 
the acoustic propagation is almost not perceptible with an acoustic ray tracing model 
as shown by Figure 5.25. A transmission loss model would be necessary to account for 
the transmission loss due to the reflection on a rough sea surface. 

3. Conclusions for June 

In June, and for the reminder of the summer, the sound propagation for a 
deep transmitter is characterized by a sound-channel type of propagation and 1s not 
influenced by the atmospheric forcing on a short time scale of a few days. 

Beewises Or tie absence Of a mixed layer and the presence of a sharp 
thermocline due to the important effect of surface heating, the detection capabilities of 
a surface-ship sonar are poor and limited to the convergence zone with ranges between 
30 and 35 km. However, strong winds can create a shallow mixed layer in which 
shallow targets could be detected as depicted in Figure 5.26, where rays emanating 


from both shallow and deep transmitters are traced. 
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Figure 5.24 June, strong wind, source at 
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Figure 5.25 June, strong wind, source at 90 meters, t=0,48 hrs, B = 81°,82°,83°,97°,98° 99°. 
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Figure 5.26 June, strong wind, sources at 10 and 90 meters 
t=0,48 hrs, = 8°. 
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EGO NCLUSION 


On a short time scale of a few hours to three days. the atmospheric conditions 
influence the acoustic propagation only for a shallow source. The ray trace emanating 
neni a deep tramsnutter 1s almost not aflected by the atmospheric forcing. This 
conclusion 1s independent of the type of propagation, depending on the time of the 
vear (deep-sound channel, convergence zone, or half duct). 

The wind is the dominant atmospheric factor. For short time scales, the coupled 
model is not verv sensitive to the other boundary conditions (solar irradiation, and 
rainfall). This 1s a favorable finding in the sense that the wind is the easiest parameter 
to measure and to forecast at sea. However, the importance of a heavy rainfall has 
been demonstrated in the case of light wind conditions. 

The initial temperature profile plavs a determinant role and has to be provided as 
accuraiely as possible. The coupled model is not very sensitive to the initial salinity 


profile. 


A. ADVANTAGES OF THE COUPLED MODEL 

The graphic output sequence from the coupled model shows qualitatively the 
acoustic propagation and its evolution under the forcing of specified atmospheric 
conditions, and an assumed initial temperature profile. 

The simplicity of the model makes it possible to be implemented on a simple 
desktop computer with some graphics capability. 

A current XBT, digitized by using a one-meter depth increment, can be used as 


Input in the model. 


B. WEAKNESS OF THE COUPLED MODEL 

The fact that the coupled model does not account very well for the operating 
frequency is inherent to a ray tracing routine. The user of such a model will have to 
simulate a beamwidth at the source corresponding approximatively to a given 
Pequency, 

The output of the coupled model does not provide such quantitative results as 
median detection range, convergence zone range, S/N ratio at the receiver, and it does 


not account for scattering and attenuation. 
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Perfect surface reflection has been assumed whatever the sea-state. An 
infinitively-deep ocean has been assumed in order to avoid the complication of bottom 
reflection. However. the model can handle perfect bottom reflection or absorption. 
Furthermore, the earth curvature has been neglected. 

The coupled model assumes horizontal homogeneity of the water mass, which we 
considered to be realistic for the horizontal ranges considered. 

Finally, we neglected horizontal advection, which is compatible with the short 


tume scales studied, especially in the Mediterranean where the currents are weak. 


C. RECOMMENDATIONS 

Further simulations could be applied to different locations in the Atlantic, 
Pacific, or Indian Oceans. where climatology is more available. 

Real data could be used to initlahze the coupled model and the boundary 
conditions could be computed from observed atmospheric and oceanic conditions. The 
resulting outputs of the model could be compared with corresponding simulations. 

Then, one might couple the OBL model with one of the different acoustic models 
used by the LS Navy (such as RAYMODE, FACT, or PE niodels) in order toneeean 
quantitative results including MDR or CZR for usual operating frequencies and a 


variety of environmental conditions. 
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APPENDIX 
RAY TRACING SUBROUTINE USE EXAMPLE 


File giving an example of sound speed profile to be used with RAY 


Or 1482. 
SO. 1452.0 
ao 1484.1 
B20. . 1480.5 
JLSO)e 1482. 
oD. 1484.1 
00. 1436.2 
c LCDR JACQUES FOURNIOL. FRENCH NAVY. 


PARAMETER (YO=50. ,B0=90. ,DB=.5,M=7,MM=6 ,N=300, NN=5000, RANGE=30. ) 
DIMENSION BET(M),CC(0:MM) , YY(0:MM) 
DIMENSION C(0:N),G(-1:N),¥C(0:N) 
DIMENSION Y(O:NN) ,Z(0:NN) 
READ(05,10)(¥Y(I),CC(I) , I=0,MM) 
10 FORMAT(F7.2,2X,F7.2) 
CALL RAY(YO,BET,BO,DB,M,CC,YY,MM,C,G,YC,N,Y,Z,NN, RANGE) 
STOP 
END 


SUB OUPEVESUSING TRE MMtl POINTS OF A STRAIGTH LINE SEGMENT 
GermiIjU@USsPROFILE < CC VS YY ) TO TRACE M ACOUSTIC RAYS ISSUED 
Be@iien SOURCE AT A DEPTH YO. 
Poe Uli ERWATERSACOUSTICS=A LINEAR SYSTEMS THEORY APPROACH" 
Bro OMerqpme se 25/ezeon ACADEMIC PRESS, ORLANDO, FLORIDA. 1985. 
Pope SURFACE AND BOTTOM REFLECTIONS. FLAT BOTTOM. 

BEtage: STARTING ANGLE OF RAY, REF VERTICAL, IN RAD. 


Gy Oy MOO) O-O) OO). OO OO. OO a Oe a FS 


M > NB OF RAYS SENT: 

BET eke Oe ae Tt hte eNGLES IN DEGREES AT THE SOURCE. 
BO Te ONTTEALSANGEE OF THE UPPERVRAY OF THE BUNDLE. 
DB s JNCREHENT OF INITIAL ANGLES IN BET. 

YO noCURCHE DERI IN SNELER.: 

Mie. toeOr Fis JN PROVIDED SSP. 

6 : PROVIDED SSP IN M/SEC. 

pind we Oh wie tis CF wi aiSeSsSr IN METERS. 

NN 7 2nSER OF RANGE. 

sd : DEPTH IN METERS AMD 

Z mOR mOloltanNeGs TRAVELLED BY A RAY IN KH. 

| oe ee ceo OF DEPTH IN MeIER. 

c (ego e talib INCRSIENT. 


2] 


G : GRADTEND TOs scr. 

ve ; ARRAY OF DEPTE SOR SLOT ice ser. 
RANGE ; MAX. RANGE OF THE PLOTIN Ki. 

K : DEPTH SIL eE. 


N MUST BE THE INTEGER VALUE OF THE MAX. DEPTH YY(MM) OF THE 
INT TISE Sske 


M1 O-G © O70 


SUBROUTINE RAY(YO,BET,BO,DB,M,CC, YY MM 3c |G Ye li 2 NN etic 


DIMENSION BET(M),CC(0:MM),YY(0:MM) 
DIMENSION C(0:N),G(-1:N),¥C(O:N) 
DIMENSION Y(O:NN) ,Z(0:NN) 

DATA PI/3.141592654/ 

¥(0)=YY(0) 

YC(9)=YY(0) 

IMIN=0 

ClG cee 


e COMPUTE THE SSP C AND ITS GRADIENT G WITH 2 METER DEPTR INGle =e 
DO 10> f=17tii 

DYY=YY(1I)-YY(I-1) 
JMAX=IFIX(YY(I)) 
GG=(Ce(1)-=ce(t-1)) bry 
DO 30 J=JMIN, JMAK-1 
G(J)=GG 
C(J+1)=C(J)+GG 

30 YC(J+1)=FLOAT(J+1) 
JMIN=JMAZ 

8, CONTINUE 


C CREATE GRADIENT AT SFC AND BOTTOM TO HANDLE REFLECTION 
G(-1)=G(0) 
G(N)=G(N-1) 


€ DO 03 J=1,M 
C03 BET(J)=90.-FLOAT(J)*0.5 
Do oe o= ten 
03 BET(J)=BO-FLOAT(J)*DELB 


C RESEARCH OF CMIN AND CMAX FOR PLOTTING 

CMIN=1550. 
CMAX=1450. 
DOU4On =o! 
CMIN=AMIN1(CMIN,C(I)) 

40 CMAX=AMAX1(CMAX,C(I)) 
CMIN=AINT(CMIN) 
CMAX=AINT(CMAX)+1. 


C27 Gig 


(2) 


PLOT THE SSP AND PREPARE PLOT FOR RAYS 
CALL COMPRS 
ChEl Cx 1 (2007) 
CALL TEK618 
CALL SHERPA('THESEOUT','8',1) 
CALL HWROT('AUTO') 
CALL BLOWUP(0.745) 
CALL PAGE(11.,5.5) 
CALL PHYSOR(1.2,1.) 
CALL ARHAZI ies a) 
CALL THKFRM(.01) 
CALL XNAME('SSP (M/S)$!,-100) 
CALL YNAME('DEPTH (M)$!',100) 
CALL CROSS 
CALL YAXANG(0.) 
CALL YTICKS(5) 
CALL XTICKS(IFIX(CMAX-CMIN) ) 
CALL GRAF (CMIN, CMAX-CMIN, CMAX, FLOAT(N),-50.,0.) 
CALL CURVE(C,YC,N+1,0) 
CALL FRAME 
CALL ENDGR(0) 
CALL PHYSOR(2.7,1.) 
CALL AREA2D(7.3,3.) 
CALL XNAME('RANGE (KM)$',-100) 
CALL YNAME('S!,100) 
CALL YNONUM 
CALL CROSS 
CALL XTICKS(5) 
CALL GRAF(0.,5.,RANGE, FLOAT(N),-50.,0.) 


DO 05 J=1,M 

Y(0)=YO 

2/0 

K=IFIX(YO) 
BETAO=BET(J)/180.*PI 
C0=C(K) 
B=SIN(BETAO)/CO 


DO 20 I=1,NN 
NNN=I 


eee werawoec® REFLECTION 
Tucano. 0) BETAO=PI-BETAO 


e CHECK IF BOLIOM REPEE CT Teh 
IF (K.EQ.N) BETAOQ=PI-BETAO 


IF ((G(K).GT.0).AND.(G(K-1).GT.0.)) THEN 
IF (BETAO.GT.PI/2.) THEN 
K=K-1 
BETA=PI-ASIN(B*C(K) ) 
Y(I)=FLOAT(K) 
Z(1I)=Z(I-1)+(COS (BETAQ)-COS (BETA) )/ (B*G(K) ) 
BETAO=BETA 
ELSE 
BETAC=ASIN(C(K)/C(K+1) ) 
IF (BETAO.LT.BETAC) THEN 
K=K+1 
BETA=ASIN(B*C(K) ) 
Y(I)=FLOAT(K) 
Z(I)=Z(I-1)+ (COS (BETAO) -COS (BETA) ) / (B*G(K-1)) 
BETAO=BETA 
ELSE 
Y(I)=FLOAT(K) 
Z(I)=Z(I-1)+2.*COS (BETAO) / (B*G(K) ) 
BETAO=PI-BETAO 
ENDIF 
ENDIF 
ELSE 
IF ( (GOK). LT.O.). AND. (GQK-1) 22120) aia 
TE (BETAOeET.PL/2.) TREN 
K=K+1 
EETA=ASIN(B*C(K) ) 
Y(I)=FLOAT(K) 
Z(I)=Z(I-1)+(COS (BETAO)-COS (BETA) ) / (B*G(K-1) ) 
BETAO=BETA 
ELSE 
BETAC=PI-ASIN(C(K)/C(K-1) ) 
IF (BETAQ.GT.BETAC) THEN 
K=K-1 
BETA=PI-ASIN(B*C(K) ) 
Y(I)=FLOAT (K) 
Z(1)=Z(1I-1)+(COS(BETAO) -COS (BETA) )/(B*G(K) ) 
BETAO=BETA 
ELSE 
Y(I)=FLOAT(K) 
Z(I)=Z(1-1)+2.*COS(BETAO) / (B*G(K-1) ) 
BETAOQ=PI-BETAO 
ENDIF 
ENDIF 
ELSE 
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IF((G(K).LT.0.).AND.(G(K-1).GT.0.)) THEN 
IF (BETAOQ.GT.PI/2.) THEN 
K=K-1 
BETA=PI-ASIN(B*C(K)) 
¥(I)=FLOAT(K) 
Z(I)=Z(I-1)+(COS (BETAO)-COS (BETA) )/(B*G(K) ) 
BETAO=BETA 
ELSE 
K=K+1 
BETA=ASIN(B*C(K)) 
¥(1I)=FLOAT(K) 
Z(1)=Z(1-1)+(COS(BETAO)-COS (BETA) )/(B*G(K-1)) 
BETAO=BETA 
ENDIF 
ELSE 
ie CCK Gin) s AND. (G(K-1).LT.0.)) THEN 
IF (BETAO.GT.PI/2.) THEN 
BETAC=PI-ASIN(C(K)/C(K-1)) 
IF (BETAO.GT.BETAC) THEN 
K=K-1 
BETA=PI-ASIN(B*C(K) ) 
¥(1)=FLOAT(K) 
Z(I)=Z(I-1)+(COS(BETAO)-COS (BETA) )/(B*G(K) ) 
BETAO=BETA 
ELSE 
Y(I)=FLOAT(K) 
Z(I)=Z(I-1)+2.*COS (BETAO)/(B*G(K-1)) 
BETAO=PI-BETAO 
ENDIF 
ELSE 
BETAC=ASIN(C(K)/C(K+1)) 
IF (BETAO.LT.BETAC) THEN 
K=K+1 
BETA=ASIN(B*C(K)) 
¥(I)=FLOAT(K) 
Z(I)=Z(I-1)+(COS(BETAO)-COS (BETA) )/(B*G(K-1)) 
BETAO=BETA 
ELSE 
Y(I)=FLOAT(K) 
Z(I)=Z(I-1)+2.*COS (BETAO)/(B*G(K) ) 
BETAO=PI-BETAO 
ENDIF 
ENDIF 
ENDIF 
ENDIF 
ENDIF 
ENDIF 


IF(Z(1I).GT.RANGE*1000.) GO TO 21 
CONTINUE 


CONTINGE 


CONVERSION M TO KM 
DO 50 I=0,NNN 
Z(1)=Z(1I)*0.001 


PLOT THE RAYS 
CALL CURVE(Z,Y,NNN+1,0) 


CONTINUE 


CALL FRAME 

CALL HEADIN( 'ACOUSTIC RAY SIRACING. = 0C7 2a) 
CALL ENDPL(0) 

CALL DONEPL 

SlCr 

END 
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